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Abstract. We compute the continuum limit of the spectra for the XX-model with arbitrary 
complex boundary fields. In the case of hermitian boundary terms one obtains the parti- 
tion functions of the free compactified boson field on a cylinder with Neumann-Neumann, 
Dirichlet-Neumann or Dirichlet-Dirichlet boundary conditions. This applies also for certain 
non-hermitian boundaries. For special cases we also compute the free surface energy. For 
certain non-hermitian boundary terms the results are more complex. Here one obtains loga- 
rithmic corrections to the free surface energy. The asymmetric version of the XX-model with 
boundaries (this includes the Dzyaloshinsky-Moriya interaction) is also discussed. 



1. Introduction 

In the first paper of this series [1] we have shown how to diagonalize the Hamiltonian 
of the XX-model with boundaries 
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h = - K +(T 7+i + a J a t+i] + 7f i a - a i + a + a t + a * a i + + + p*<ri] (!) 
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by introducing an auxiliary Hamiltonian 
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+ 7i + a+(J o (7 i + + a ^ + P+ a i a L+i + P-^i+i + Pz° z l] ( 2 ) 

which in turn may be diagonalized in terms of free fermions. The parameters a±, j3±, a z 
and (3 Z are arbitrary complex numbers. Note that H\ ong commutes with cJq and cx£ +1 . 
Hence the spectrum of Hy ong decomposes into four sectors (+, +), (+, — ), (— , — ), (— , +) 
corresponding to the eigenvalues ±1 of (7q and cr£ +1 . The spectrum of H is obtained by 
projecting onto the (+, +)-sector. 

While in [1] our focus was on the diagonalization of the finite chain, here we are 
going to examine the asymptotic behaviour of the energy levels for large values of L. 
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In the case of conformal invariance this information is usually encoded in the partition 
function, which we define by 

Z = lim trzi {H ~ e ° oL ~ f ° o) (3) 

L^oo 

where z < 1 and and denote the free bulk and the free surface energy respectively. 
£ may be considered as a normalization constant here. We will also consider the partition 
function for H\ ong . 

In the case of periodic boundary conditions the continuum limit of the XXZ-chain 
can be described by a free boson field $ (see e.g. [2]), with action 

S = i J d Xl dx 2 [{d^f + (<9 2 $) 2 ] (4) 

being compactified on a circle of radius r, i.e. we identify $ with $ + 2nr (we follow 
the notation of [3]). The radius of compactification r is related to the value of the 
anisotropy in the XXZ-model. At the free fermion-point the anisotropy is zero and we 
have r = 1 / \f^ix. 

Consider the boson field on a cylinder of length L and circumference N imposing 
periodic boundary conditions in the direction of the cylinders length. Quantizing the 
field with time in iV-direction yields the partition function [3] 

Aperiodic = ^T^ \qA^ ) 5 

where q = e - 27rN / L an d 

oo 

fl(q)=q*l[(l-q n ). (6) 

n=l 

This expression coincides with the partition function for the periodic XXZ-chain, which 
has been studied extensively in [4, 5, 6]. The partition function for H\ ong in the periodic 
case is just 4Z pcriodic . 

Suppose now the case of open boundaries at both ends of the cylinder. There are 
two different boundary conditions preserving conformal invariance, i.e. the Dirichlet 
and the von Neumann (see e.g. [3, 7] ). This yields three different partition functions 
corresponding to the various possible combinations of the two boundary conditions. 
The respective partition functions can be found in [3, 7]. Imposing Dirichlet boundary 
conditions at both ends of the cylinder results in 

iWA) = j-£g^™+ A ) 2 (7) 

where q = e~ wN ^ L and A = ($ ~ ^l)/Vtt- By $ an d $l we denote the values of the 
boson field at the boundaries. This type of partition function has also been obtained 
for the open XXZ-chain with diagonal, hermitian boundary fields [8]. 
The Neumann-Neumann partition function is given by 

Z^) = ^Y,^ m ' {2 ^ r)+N2? ( g ) 
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where A = (<3> — ^l)/V^ and $ , $ L denote the values of the dual field <3> at the 
boundaries [3, 7]. The Dirichlet-Neumann boundary condition yields [3, 7] 



In this case no free parameter appears. 

Which boundary conditions of the Hamiltonian H lead to the partition functions 
(8) and (9) is yet not known. We will close this gap. We are going to see that, as long 
as the Hamiltonian H is hermitian, the partition function of the chain corresponds to 
one of the three boson partition functions for open boundaries. However, this may also 
be the case for certain non-hermitian boundary terms. See section 2 for a survey over 
the boundary conditions which we will consider. 

In [9] it has been argued that non-diagonal boundary terms correspond to von 
Neumann boundary conditions for the boson field. Our results will show that this 
assumption is correct in the hermitian case, but not necessarily for non-hermitian 
boundaries. 

In this paper we also study the asymmetric XX-chain with boundaries (see equation 
(25)). This includes the quantum chain with Dzyaloshinsky-Moriya interactions, which 
was already studied in [10] for periodic boundary conditions. The partition function for 
real values of p and q in (25) and periodic boundary conditions was obtained in [13]. 
The special case p — 1 and q = with boundaries has already been studied in [11]. 

This paper is organized as follows: We will start with a summary of our results for 
the symmetric chain with hermitian and non-hermitian boundaries in section 2. The 
asymmetric XX-chain will be discussed in section 3. We will give our conclusions in 
section 4. The appendix is dedicated to the computation of the fermion energies which 
yield the energy gaps of the chain. 

2. Summary of results 

The structure of our results is most simply encoded in terms of new parameters 



F, C, G, K, J, where 

F = Aa-a + (3 + /3- 
C = a 2 _Pl + a 2 + (3 2 _ 

G = 2a_a+(l + 2f3 2 z ) + 2/3_/3+(l + 2a 2 z ) 
K = 2(a_a+ + (3_(3 + ) + (1 + 2a 2 z )(l + 2{3 2 z ) 

J = 8a z p z - Aa 2 z p 2 + 2(a_a+ + (3_(3 + ) + 2a 2 + 2(3 2 z - 1. (10) 



The six different situations we considered are given in table 1. This table also indicates 
which kind of partition function we obtained for H. For hermitian boundaries the cases 
in this table cover the whole parameter space. This is not the case for non-hermitian 
boundary terms. 




(9) 
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Table 1. Conditions on the parameters F, C, G, K, J being considered in this paper. 
An asterisk indicates that the corresponding value may be arbitrary. The partition 
functions Zdd,Znm and Zdn are defined in (7), (8) and (9). For the last two cases 
we obtained anomalous behaviour of the energies (see section 2.5). We use h. and n-h. 
as a shortcut for hcrmitian and non-hermitian. 

F C G K J 



^ * * * * Z NN h. and n-h. 

^ * * Z DN h. and n-h. 

^ * Z DD h. and n-h. 

2 x Z DD n-h. 

^0^0 * * anom. n-h. 

^0 anom. n-h. 



For the cases in table 1 we computed analytically the energy gaps in the leading 
order as shown in appendix A. The exact ground-state energies for certain cases being 
listed in table 2 have already been obtained in [1]. 

However, the information of the energy gaps is already sufficient to obtain the 
partition functions for H\ ong upto a factor z ~^ +hmin , where h min denotes the lowest 
highest weight appearing in the representation of the Virasoro algebra. This information 
has to be extracted from the ground-state energies [12] 

E = e oo L + f x -j(-^-h mi ^\+... (11) 

The value of the free bulk-energy is already known from the periodic chain, i.e. 600 = —^ 
[14]. The values of the free surface energies and of h m [ n have been obtained for the 
cases given in table 2 expanding the exact expressions for the finite chain given in [1]. 

The results confirm the partition functions we will give for ifi ong in the following. 
For the cases not being listed in table 2 we checked our expressions numerically. The 
partition functions for H have been obtained afterwards by projecting onto the (+, +)- 
sector (see [1] for details on the projection mechanism). We are now going to discuss 
the cases given in table 1 in detail. 

2.1. Neumann-Neumann boundary conditions 

The condition F^O simply implies that all non-diagonal boundary terms are present. 
The diagonal terms may be arbitrary. It is quite instructive to introduce a new 
parametrisation of the boundary parameters a± and f3± as follows: 

a+ = Rj^ a_ = R a e^ f3 + = iV (x+v) (3- = Rpe~' m{x+ip) (12) 

where R a , Rp G 1R + , x, y? G R, — 1 < x < 1- The value of A which enters the partition 
function (8) is then given by 

AH X forLodd (13) 
1 x + 1 for L even 



Table 2. The values of the parameters F, C, G, K, J for the cases where the ground-state energy of i?i ong has been obtained in [1]. 
The table shows the free surface energy /oo and the lowest highest weight /i m ; n (cf. (11)). The expressions for are also valid for H. 
C, denotes a free parameter, which has to be chosen such that < Re£ < 2n for case 12 and — 7r < Rc( < tt otherwise . 
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The partition function for H\ ong is given by 

Z long = 2Z NN ( X ) + 2Z NN ( X + 1) (14) 

for even and odd values of L. Note that the parametrization (12) does not work for 
non-hermitian boundaries. In this case we may define the parameter x via 

«-5 k (S) (15) 

where the logarithm is taken in a way such that < |Re(%)| < \ for Re(o; + /9_ + a_/5 + ) > 
and 1 > |Re(x)| > \ for Re(a +( 5_ + a_f3 + ) < 0. This rule is consistent with the 
parametrisation (12), which implies Re(a + /3_ + = 2R a Rp cos(x7r). However, we 

cannot perform the projection onto the (+, +)-sector in general if the boundaries are 
non-hermitian. In this case our result for the partition function of H is restricted to the 
cases listed in table 2 which satisfy the conditions a z = j3 z = or a + = f3 + = f3- 
[1]. The result for H\ ong is valid in general. 

2.2. Dirichlet- Neumann boundary conditions 

The conditions for the Dirichlet-Neumann case yield 

/3± = a± ^ (3 Z ^ ±-^= or /3± ^ a± = a z ^ ±-J=. (16) 

v 2 v 2 

The value of a z respectively of j3 z may be arbitrary. The partition function for H\ ong is 

just Z\ ong = 4Zdn for L even and odd. The case of one non-diagonal boundary at the 

end of a semi-infinite XX-chain has already been considered in [15]. The author has 

shown, that the model can be decoupled into two Ising models with different boundary 

conditions. One of them being subject to a free boundary condition at the end, the other 

one being subject to the fixed boundary condition. The partition functions for the Ising 

model with free and mixed boundary conditions are well known [12, 16]. Taking the 

product of them also results in (9). 

2.3. Dirichlet- Dirichlet boundary conditions 

The Dirichlet-Dirichlet partition function is obtained for two different types of boundary 
terms. The first type is given by 

&^±-J= ^ z ^±^= (3_ = a_ = 0. (17) 

where instead of the last equation we might have also chosen f3 + = a + = 0. However, 
for this type of boundaries the non-diagonal boundary terms have no influence on the 
energies even on the finite chain. The second type is given by 

/3 ± = a± ^ /3 Z = ±-^= or f3± ^ a± = a z = ±-^= (18) 

v 2 v 2 
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where the values of a z respectively (3 Z may be chosen arbitrarily. Note that for this 
second type the boundary terms are non-hermitian. For both types of boundaries the 
value of A which enters (7) is given by 



A = — arccos 
2vr 



J 



(19) 



The partition function for H\ ong is Z\ ong = 4Zdd(A) for even and odd L. 



2.4. F = C = G = K = J = 



For the case F = C = G 



a. 



±- 



: K 



J = the boundary parameters are nearly fixed, i.e. 
i 



V2 



(20) 



V2 

and a + = f3 + = or a_ = /3- = 0. This case is special. The partition function we 
obtained does not fit into the picture suggested by the analogy to the free boson. Here 
an additional zero mode appears. The value of A is given by A = for odd L and by 
A = 1/2 for even L. The partition function for ifi ong is Z\ ong = 8Z DD (A) for L even 
and L odd. 



2.5. Anomalous behaviour 

For the last two cases in table 1 we obtained logarithmic corrections to the free surface 
energy. Note that there is no case in table 2 which corresponds to the boundaries in 
question. However, we computed the ground-state energy E of H\ ong numerically for 
one example of boundaries which lead to the anomalous behaviour. Our computations 
suggest 



E ~ -- + foo + \ (pi In L + p 2 ) 

7T L 



(21) 



where pi,p2 are complex numbers. Figure 1 and figure 2 show the real respectively 
the imaginary part of the Casimir amplitude. The computation has been done for 
a z = 1/V2, f3 z = i/y/2,a± — @± — which corresponds to the last case in table 1. In 
both cases the fermionic energies scale as 



2A ~ - 

Lj 



kn ± - (arg A 



ilnL) 



k = 0,1,2,.... 



where the value of A is given by 

A = -^ 
G 

for F = 0, C, G 7^ respectively by 

A = (-ir 2J 



(22) 



(23) 



(24) 



1 - 4o2# 

for F = C = G = K = 0, J 0. Our result for the last case makes only sense for 
(a z (3 z ) 2 7^ 1/4. The factor 2 in (22) is just a remnant from our notation in [1]. Note the 
integer spacing of the energy gaps (22). 



Rc[L(E a + L/n + 0.25)} 

2 




1000 



Figure 1 . Real part of the Casimir amplitude (a z = 1 / \/2, f3 z = i/ \/2, ct± = Pi 



lm(LE ) 

0.75 
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Figure 2. Imaginary part of the Casimir amplitude (a z 
/3±=0). 



- i/V2,a± 
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3. The asymmetric XX-chain 

We also considered the asymmetric XX-chain which is defined by the Hamiltonian 

+ -j= [o^f + + a^ 2 + ff_a L + + p z a z L ] . (25) 

Without loss of generality we restrict ourselves to ^fpq = \. Under this condition H & 
can be mapped on the symmetric chain (1) with boundary parameters 

a_=Q^a'_ a + = Q^a' + (3_ = (3'_ f3 + = Q^r (3' + (26) 

by a similarity transformation, where Q = The diagonal terms remain unchanged. 
We will separate two cases in the following. First we will consider hermitian bulk terms, 
i.e. 

p = e inK /2, q = e~ i7r 72. (27) 

This corresponds to Dzyaloshinsky-Moriya type interactions in the bulk. Thereafter we 
will turn to non-hermitian bulk terms with p, q in R. 

3.1. Dzyaloshinsky-Moriya interactions 

We restricted ourselves to hermitian boundaries. In this case the value of k has only an 
effect onto the spectra if non-diagonal boundary terms are present at both ends of the 
chain. We introduce the parametrization 

a' + = R'J^' ol_ = R' a e~' mip ' (3' + = R'^'^ (3'_ = i^e"^'^ .(28) 

According to (26) and (12) the mapping onto symmetric bulk terms yields 

R' a = R a R'p = R P <p = v' + k^^ X = X +<L-l)mod2. (29) 

Remember that in the case of length independent boundary terms, the value of \ 
contains the full information about the partition function. Observe now, that it is 
possible to choose values of k such that the value of x m (29) becomes independent of 
L as long as one considers only certain sequences of lattice lengths. This is exactly the 
case if k/2 is rational, i.e. k/2 = m/n,m e Z, n G N + . In this case x is independent of 
L for L = In + r, < r < n, I G N, i.e. 

n 

The results obtained for the symmetric case can be adopted immediately. This is not 
possible for irrational values of k/2. Note that this kind of commensurability and 
incommensurability has also been observed for the periodic chain with this type of 
interaction [10]. 
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3.2. p,qeR 

Without loss of generality, we are going to restrict ourselves to p > q which implies 
Q < 1. We have to distinguish two different situations. If a'_(3' + equals zero we can 
adopt the results for the symmetric case with F = C = (see appendix A for details). 
One has just to exchange a± and /3± by a'± and If otherwise a'_j3' + ^ the situation 
changes. The energy gaps suggest the partition function for the long chain (2) with 
boundary terms given by (26), i.e. 



where 



and 



r>(z) 



r = (1 - Q 2 (l - 2a 2 z - 2a'_a' + ) - 2a 2 Q 4 ) 

x (1 - Q 2 (l - 2(3 2 z - ~ 2(3 2 z Q')/(a' 2 (3' 2 Q\l + Q 2 ) 2 ). (33) 

Note that the value of A y is purely imaginary. Hence the length dependent term given 
in the partition function is a phase. Such a term already appeared in the toroidal 
partition function for the asymmetric model with periodic boundary conditions [13]. 
Note also that the result (31) simplifies to the partition function we obtained in the 
Neumann-Neumann case for i?i on g if one sets Q — 1 (see (14)). 

In [1] we computed the exact ground-state energy on the finite chain for boundaries 
defined by 

a z = (3 Z = a' + a'_ = (3' + (3'_ = 1. (34) 
If we introduce the parameter x' via P' + = e mx ' a' +) then (33) simplifies to 

T=(e mx 'Q 2 y 2 . (35) 
Expanding the exact expression obtained in [1] leads to 

q + q- 1 Q + Q- 1 + (Q - Q-^inr (Q + Q-y /i (inrn 

~ 2tt 2tt 2L V 24 8tt J ' { ' 

It is only for this case that we can perform the projection onto the (+, +)-sector. We 
obtain 

Z = J_Y z 2( m+X >/2~nnQMl + 2mLA y for L ,3- 

Z=^— Y z 2(m +X >/2-nnQMZ + 2 m LA y f Qr even £ ( 38 ) 
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Note that our results for i^i on g only apply if (33) is different from zero. If the nominator 
of r vanishes we obtain logarithmic terms in the asymptotic behaviour of the energy gaps 
similar the ones obtained for the symmetric chain. We obtained the fermion energies 



2A ~ 
where 



Q + Q 
2L 



-i 



kn ± - 
2 



arg A — i In L 



±i 



.Q-Q- 



2a'l(3'lQ\l + Q 



2\2 



(39) 



(40) 



2AQ 2 + AQ\B + E 2 ) - QQ 6 (D + 2E 2 ) + 8E 2 Q 8 ' 

The values of A,B,D,E are given in (A. 3). We are not going to consider the case, 
where the denominator in (40) vanishes. 



4. Conclusions 



In this paper we considered the spectra of the XX-model with boundary fields given by 
the Hamiltonian H in (1). In order to obtain our results we also studied the Hamiltonian 
.f/iong defined by (2). Furthermore we considered the asymmetric XX-chain (25). Here 
we separated two cases. First we considered the case where the values of p and q in 
(25) are given by p = e 17m /2 and q = e~ 17rK /2 (this corresponds to Dzyaloshinsky-Moriya 
interactions). Second we considered the case where p and q are real numbers. 

For periodic boundary conditions the partition function for H is given by the 
partition function of the free boson with periodic boundary conditions (5) [4, 5, 6]. The 
partition function for H\ on% with periodic boundary conditions is just the expression 
(5) multiplied by 4. The spectra of the asymmetric chain with Dzyaloshinsky-Moriya 
interactions and periodic boundary conditions have been studied in [10], whereas the 
partition function for real values of p and q and periodic boundary conditions was given 
in [13]. 

The results we obtained for H and H\ ong can be encoded in terms of the parameters 
given in (10). For reasons discussed in the text we obtained the energy gaps only for the 
cases given in table I. Note that as long as we restrict ourselves to hermitian boundary 
terms we studied the most general case. This is not true for non-hermitian boundaries. 

For H with hermitian boundaries we obtained the partition functions corresponding 
to one of the three boson partition functions (7), (8) and (9). However, we found these 
partition functions also for certain non-hermitian boundary terms. 

We obtained the Neumann-Neumann partition function (8) if all non-diagonal terms 
are present, where the value of A in (8) is given by (13). The Dirichlet-Neumann 
partition function (9) is obtained for boundary terms which satisfy (16). We found the 
Dirichlet-Dirichlet partition function (7) for two types of boundary terms given by (18) 
respectively (17). The value of A in (7) is given by (19) in both cases. Furthermore, for 
the case of non-hermitian boundaries we found a case which is special. For boundary 
terms given by (20) the partition function is the Dirichlet-Dirichlet partition function 
(7) multiplied by 2, where the value of A is or 1/2 for odd respectively even values of 
the lattice length. 
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The partition functions for H\ ong have also been considered. If all non-diagonal 
terms are present it is given by (14). For the other cases it is just the partition function 
of H multiplied by 4. We also computed the values of the free surface energies of H\ ong 
and H for certain boundary terms. They are given in table 2. In this table we have also 
given the values of the lowest highest weights which appear in the spectra of H\ ong for 
these boundary terms. 

For the last two cases in table 1 we obtained logarithmic corrections to the free 
surface energy (equation (21)). Here we found also logarithmic terms in the asymptotic 
behaviour of the energy gaps (equation (22)). 

In the case of Dzyaloshinsky-Moriya interactions, we restricted ourselves to 
hermitian chains. The value of the phase k (see (27)) has only an effect on the spectra 
if non-diagonal boundary terms are present at both ends of the chain. In this case we 
obtain the Neumann- Neumann partition function (8) if k/2 is a rational number. The 
value of x m the definition of A in (13) has just to be exchanged by the expression in 
(30). 

The partition function for the asymmetric chain (25) with real values of p and q has 
only been obtained for one special type of boundaries (see equation (34)). It is given 
by (37) for L odd respectively by (38) for L even. The asymptotic behaviour of the 
ground-state energy for this case is given in (36). 
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Appendix A. Determination of the energy gaps 

In [1] we have seen that the spectrum of Hi ong is given in terms of free fermions. The 
fermionic energies 2A n are given by 

2A„ = l^ + rr; 1 ). (A.l) 

where the the roots of the polynomial 

p{x>) = L4L+8 + 1 _ + ^ + (£? + + ^ 

+ (D + 2E 2 )(x iL+2 + x 6 ) + E 2 (x 4L + x 8 ) - 2E(x 2L+8 + x 2L ) 
+ {{A - B - D - l)/2 - (-1) L C - 2E 2 ) (x 2L+6 + x 2L+2 ) 

+ (A - B - D - 1 + 2(-l) L C + AE- AE 2 ) x 2L+4 } I (x 2 - l) 2 (A.2) 
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where the coefficients appearing in (A. 2) are functions of the boundary parameters: 

A = 2(a_a+ + + a 2 z + (3 2 z - 1) C = a 2 _{3 2 + + a\(3 2 _ 

B = (2a_a + - l)(2/3_/3+ - 1) + 4(3 2 z {a_a + - 1) + 4a 2 z {(3 + (3_ - 1) 

D = (3 2 (4a.a + - 2) + a 2 z (4{3 + {3_ - 2) E = 2a z (3 z . (A.3) 

Since there appear 4 zeros for each fermion, namely x ni x~ 1 ,— x n and —a;" 1 , the 
polynomial yields L + 1 fermionic energies (see [1] for details). In addition to these 
fermions there always exists a fermion with energy 2A = 0, which we named 'spurious' 
zero mode. All possible combinations of these L + 2 fermions build up the spectrum 
of .f/iong- The spectrum of H is then obtained by excluding the 'spurious' zero mode 
from the set of fermion energies and then taking the sector with an even or with an odd 
number of fermions being excited with respect to the vacuum. How to decide whether 
one has to pick up the even or the odd sector has been discussed in detail in [1]. 

Since in this paper we are interested into the large L behaviour of the low lying 
energy levels, we look for the asymptotics of the zeros of p(x 2 ) in the vicinity of the 
point x — ±i using the ansatz x = e 1 ^ -1 "^ , where we assume to be a constant in leading 
order. For the first four cases given in table 1 we obtained the following equations: 

o F ^ 

2C 



o F = C 



o F = C 



cos(20) + 

0sin(20) : 
G = 0,K ^ 

2 







\L+1 



J 



K 







(A.4) 



(A.5) 



(A.6) 



o F = C = G 



cos(20) + (-1) J 
K = J = 

4 exp(2i0) = (-1) L (A.7) 

Solving these equations for yields the energy gaps (cf. (A.l)) in leading order. We 
obtained the following expressions: 

o F ^ 



2A~ \ 



2n-l 1 

7r ± -arccos 

2 2 



2C 

T 



1 < n 



o F = C = 0,G ^0 
2A r 

o F = C = G 



1 n 



L2 
0,K ^ 



-7r < n 



2A ~ - 

1j 



2n-l ,1 

7i ± -arccos 

2 2 



J 



1) L+1 — 
' K 



1 < n 



(A.8) 



(A.9) 



(A.10) 



In addition to these modes there appears an additional zero mode. 
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o 



F= C=G=K=J=0 



2A ~ ymr 1 < n for L even (A. 11) 

1 2n — 1 

2A ~ - — - — 7r 1 < n for L odd (A.12) 

Furthermore we found 3 additional zero modes for L even and 2 additional zero 
modes for L odd. 

For the last two cases in table 1 the asymptotic zeros of the polynomial can only be 
found using the ansatz x = e 1 ? -1 ^ if we assume the imaginary part of to diverge as 
L goes to infinity. Hence we generalize our ansatz to x — e 1 ^ -1- ^, where (j>(L) is a 
complex function. Furthermore we assume lim^oo = and lim^oo e~ 1( ^ L ) = 0. 
The second assumption will be explained shortly. In both cases our ansatz leads to the 
solution of an equation of the form 



\ L J 



0. (A.13) 



L 

From this expression one can see that our second assumption is indeed necessary to 
solve this equation, since otherwise the first term would be finite for all values of L, 
whereas the second term vanishes as L goes to infinity. Neglecting the terms of order 
4>(L)/L this equation is solved in terms of the so called Lambert W function C which is 
defined by the property C(x)e c ^ = x. We obtain 

<j>(L) = -i£(AL)/2 where A = ~. (A.14) 

b 

The asymptotic behaviour of C is well known [17], i.e. 

C(L) ~ 2ikn + InL - ln(hiL + 2ink) + ... (A.15) 

Note that this expression is in accordance with our assumptions we made concerning 
the asymptotic behaviour of <f>(L). 

In this paper we also considered the asymmetric chain, which is similar to the 
symmetric chain with length dependent boundary parameters (cf. (26)). This length 
dependence enters the polynomial (A. 2) only via the coefficient C, which becomes 

C = a' 2 _f3'lQ 2 " 2L + a'l!3 a Q 2L - 2 . (A.16) 

For Q < 1 (the case considered in section 3.2) the first term on the RHS of (A.16) 
diverges exponentially whereas the second term vanishes exponentially as a function 
of the lattice length L. Hence we have to distinguish two different situations. If a'_(3' + 
equals zero we may find the asymptotic zeros of the polynomial as for the symmetric case 
with F = C = 0. If otherwise a'_(3' + ^ then C diverges as L is increased. This has to 
be compensated by modifying our ansatz to x — e 1 ^ _1 i +lnQ . This works as long as the 
nominator of T in (33) is different from zero. Otherwise we have to modify our ansatz 
another time to x — e 1 ^ -1- ^ hln< 3 ) where lim^oo = and lim^oo e -1 ^^^ = 0. This 
leads again to an equation of the form (A.13) which can be solved as described above. 
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